Investigating superconductor-insulator transition in thin films using drag resistance: 

Theoretical analysis of a proposed experiment 
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The magnetically driven superconductor-insulator transition in amorphous thin films (e.g., InO, 
Ta) exhibits several mysterious phenomena, such as a putative metallic phase and a huge magnetore- 
sistance peak. Unfortunately, several conflicting categories of theories, particularly quantum-vortex 
condensation, and normal region percolation, explain key observations equally well. We propose 
a new experimental setup, an amorphous thin-film bilayer, where a drag resistance measurement 
would clarify the role quantum vortices play in the transition, and hence decisively point to the cor- 
rect picture. We provide a thorough analysis of the device, which shows that the vortex paradigm 
gives rise to a drag with an opposite sign and orders of magnitude larger than the drag measured if 
competing paradigms apply. 



The superconducting state and the metallic Fermi- 
liquid form the very basis of our understanding of corre- 
lated electron systems. Nevertheless, the transition be- 
tween these two phases in disordered films is shrouded in 
mystery. Experiments probing this transition in amor- 
phous thin films such as Ta, MoGe, InO, TiN, etc., used a 
perpendicular magnetic field and disorder (tuned through 
film thickness) to destroy superconductivity. But instead 
of a superconductor-metal transition, they observed in 
many cases a superconductor- insulator-transition (SIT) 1 . 
The "dirty boson" model 2 propounded the notion that 
the insulator is the mark of vortex condensation, and 
that the SIT occurs at a universal critical resistance, 
Ru = h/4e 2 . More recent experiments, however, showed 
the critical resistance to be non-universal 3 . Furthermore, 
in many field tuned experiments, a surprising metal- 
lic phase intervenes between the superconductor and 
insulator 4-6 , with a temperature- independent resistance 
below T ~ 50mK, and (at least in Ta films) a distinct 
nonlinear I — V characteristics 7 . Quite generically 5 ' 6 ' 8 , 
these films exhibit a peak in the magnetoresistance(MR) 
curve (particularly strong in InO and TiN) as in FIG. la. 

Two competing categories of theories may account for 
these phenomena. On one hand, within the quantum vor- 
tex pictures 2, ,10 ; the insulating phase implies vortex con- 
densation, the intervening metallic phase is described as 
uncondensed vortex liquid (e.g., vortex Fermi liquid), and 
the high field nonmonotonic MR indicates the appear- 
ance of a finite electronic density of states at the Fermi 
level. On the other hand, the percolation paradigm 11 ' 12 
describes the films as consisting of superconducting (SC) 
and normal puddles; at the MR peak SC puddles exhibit 
a Coulomb blockade, and the percolating normal regions 
consist of narrow conduction channels. Yet a third the- 
ory tries to account for the low field SC-metal transi- 
tion using a phase glass model 13 (see, however, Ref. 14 
which argues against these results), but does not address 
the full MR curve. Qualitatively, both paradigms above 
are consistent with MR observations, and recent tilted 
field, AC conductance, Nernst effect, and Scanning Tun- 
neling Spectroscopic measurements 15 cannot distinguish 



FIG. 1: (a) A typical magnetoresistance (MR) curve of amor- 
phous thin film superconductors. As the magnetic field B 
increases, the superconducting phase is destroyed, and a pos- 
sible metallic phase emerges. After which the system enters 
an insulating phase, where the MR reaches its peak. The re- 
sistance drops down and approaches normal state value as B 
is further increased, (b) Our proposed bilayer setup for the 
drag resistance measurement. A current bias I\ is applied in 
one layer, and a voltage V2 is measured in the other layer. 
The drag resistance Rd is defined as Rd = Vijl\. 



between them. 

Can we design an experiment that qualitatively distin- 
guishes between the two paradigms? Here we propose a 
thin film " Giaever transformer" 16 as such an experiment 
(FIG. 1(b)). The original design of a Giaever transformer 
consists of two type-II superconductors separated by an 
insulating layer in perpendicular magnetic fields. A cur- 
rent in one layer moves the vortex lattice in the entire 
junction, yielding the same DC voltage in both layers. 
Determining the drag resistance Rd — V2/I1 in a sim- 
ilar bilayer structure of two amorphous superconduct- 
ing thin films should qualitatively distinguish between 
the two paradigms (see also Refs. 17 ): Within the vortex 
paradigm, vortices in one layer drag the vortices in the 
other, but within the percolation picture, the drag re- 
sistance is solely due to interlayer "Coulomb drag", as 
studied in semiconductor heterostructures 18 . The sign 
and the magnitude of the drag within the two paradigms 
are different: vortex drag implies the same sign for the 
voltage drops in the two layers, sign(Vi) = signiVz), but 
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the Coulomb drag yields an opposite sign for V2 and V\. 
In addition, a vortex drag would be much stronger than a 
Coulomb drag, because the films' charge carrier density is 
orders of magnitude larger than the vortex density; drag 
effects are typically inversely correlated with carrier den- 
sity. Indeed, we find that two identical films as in FIG. 
2b of Rcf. 5 with 25nm center-to-center layer separation 
at 0.07K would produce a drag resistance ~ O.lmfi ac- 
cording to the vortex theory, but only ~ 10~ 12 f2 for the 
percolation theory (FIG. 2). Below we will support these 
claims by analyzing the drag in the thin film symmetric 
bilayer within a representative theoretical framework in 
the vortex 10 and percolation paradigms 12 . 
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FIG. 2: Drag resistance Rd (in Ohms) between two identi- 
cal films as in FIG. 2b of Ref. 5 (a) vs. magnetic field B, 
according to the vortex picture 10 (log scale); (b) vs. normal 
metal percentage p (corresponding to magnetic field), accord- 
ing to the percolation picture 12 . The drag resistance in (a) 
has been smoothened to avoid discontinuity at the boundary 
between the metallic and the insulating phase. Center-to- 
center layer separation a = 25nm, temperature T = 0.07K 
and 0.35K. Insets: single layer magnetoresistance (MR, log 
scale) reproduced in each theory. The parameters are tuned 
to make the MR resemble the experimental data in FIG. 2b 
of Ref. 5. In the quantum vortex picture, Rd has a peak 
at the steepest point (~ 8T) of the MR, which is due to the 
fact that Rd is proportional to the square of the slope of the 
MR in the small magnetic field side of the peak. Also, Rd 
is larger at lower temperature, because the MR curve is then 
much steeper. Carrying out the experiments at even lower 
temperatures may further enhance the vortex drag effect. In 
the percolation picture, the sign of the voltage drop of the 
passive layer is opposite to that of the driving layer, and the 
maximum magnitude value of Rd is much smaller, ~ 10 — 12 £1. 



Drag resistance, Rd, in the quantum vortex picture. 
Within the quantum vortex paradigm, the insulating 



phase has been explained as a superfluid of vortices by 
the "dirty boson" model of Rcf. 2, while the metallic 
phase is expected to be an uncondcnscd vortex liquid. 
This picture has been pursued by Ref. 10 which argues 
that vortices form a Fermi liquid for a range of mag- 
netic field, thereby explaining the metallic phase. At 
fields larger than the insulating phase value, it is claimed 
that spinons (unpaired fermions with finite Fermi-energy 
DOS) become mobile, impede vortex motions, destroy 
the insulating phase, and suppress the resistance down 
to normal metallic values. 

Without interlaycr Joscphson coupling, the vortex 
drag comes from the magnetic coupling between vortices 
in different layers which tends to align themselves ver- 
tically to minimize the magnetic energy. To calculate 
the vortex drag, we follow the vortex-boson duality for- 
malism of Ref. 2 and include the effect of physical elec- 
tromagnetic field to obtain the Lagrangian for vortices, 
which features the following form of the intralayer vortex 
interaction potential Ui and the interlayer one U e : 



Ui(q) = 
U e (q) = 



4>l<ic 



Qc 



2tt q{q 2 + 2q c q + q^(l-e- 2 i a )) 
-e-« a Ui{q)q c /(q + q c ), 



(1) 



where q c = d/(2\ 2 ) is the inverse Pearl penetration 
depth, d is film thickness, and a is the center-to-center 
interlayer distance (we verified that accounting for the 
finite thickness gives roughly the same results as sim- 
ply taking the interlayer distance to be a center-to-center 
one) . q c can be determined from the Kosterlitz-Thouless 
temperature, Tkt, of the sample; typically q~ x ~ 1cm. 
When r < l/q c , Ui(r) gives the familiar log interaction; 
for r > l/q c , Ui{r) is still logarithmic but with half of 
the magnitude 19 , in contrast to the 1/r behavior of the 
single layer case (which is Eq. (1) with a — > 00). The 
interaction between two vortices with the same vorticity 
in different layers is attractive as expected, although its 
strength is suppressed with increasing distance a and de- 
creasing q c . These forms of vortex interaction potentials, 
which can be also derived classically by solving London 
equations for two vortices, are simply due to the mag- 
netic energy and the superfluid kinetic energy of a SC 
thin film, and thus is robust against model details. 

Deep in the insulating phase, i.e., the vortex superfluid 
phase where the vortex dynamics is presumably nondis- 
sipative, we find the drag resistance using the bilayer 
supercurrent drag mechanism 20 , applied to the vortex 
condensate. Here, a vortex "supercurrent" j\ in the first 
layer, produces a vortex "supercurrent" j 2 in the second 
layer even without interlayer tunneling. To see this, we 
use the dual vortex theory of Ref. 2 applied to a bilayer, 
and neglect dual gauge field fluctuations, which are sup- 
pressed in the insulating phase. Without vortex current 
bias in either layer, we can write a Hamiltonian describ- 
ing a vortex-superfluid in each layer, with density-density 
interaction given by Eqs. (1), and use the Bogoliubov 
approach to diagonalize the Hamiltonian. With a per- 
turbation term Hi = X}a rr Wi ' describing a vortex 
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current bias with velocity V\ in the first layer, we can 
perturbatively find the new ground state and therefore 
the drag vortex current ji in the second layer. From the 
ratio of the vortex currents in the two layers, which gives 
the ratio between the voltages in the two layers, we find 
the drag resistance in terms of the single layer resistance. 
A straightforward, but lengthy calculation yields: 



Rd 
~R~ 



32 

h 



h 



128a 2 (f) V 2i:m v n% ' 



(2) 



where n v — B/<j>o, m v are the vortex density and mass. 
We expect that in this phase the dissipative response of 
vortices, if present, is irrelevant or insufficient to local- 
ize vortices, e.g., as in the Caldeira-Leggett model of a 
dissipative Josephson junction 21 . 

In the vortex picture, the intervening metallic phase 
is interpreted as a liquid of uncondensed vortices, e.g., 
in Ref. 10, and the vortices have dissipative dynamics. 
As long as vortex superfluid is absent and the vortex 
response is dissipative, one can derive the following form 
of the the drag resistance Rd using either the Boltzman 
equation or diagrammatic techniques, irrespective of the 
effective statistics of vortices 22,23 : 



Rd — 



^r^Wo qdq L MU] ^w> 

(3) 

where U = UJ [(1 + t/ M Xi)(l + ^,2X2) - ^1X2] is 
the screened interlayer interaction, xi,2 are the vortex 
density response function of each layer. Remarkably, the 
drag resistance is proportional to dRi^/dB, which equals 
da v /dn v with a v and n v being respectively the vortices' 
conductance and density. Thus Rd peaks when the MR 
attains its biggest slope. da v /dn v appears since Rd is 
related to the single layer rectification function, T, de- 
fined as j v = T(j) 2 , with <fr being the vortex potential 
field. T is generally proportional to da v /dn v (see Ref. 
23). The only model-dependent input is the density re- 
sponse functions Xi,2- As one choice of xi,2, we follow the 
vortex Fermi liquid description for the metallic phase 10 
and use the Hubbard approximation form for xi,2 con- 
sidering the short-range repulsion between vortices and 
also the low density of this vortex Fermi liquid 24 . The 
maximum of Rd we obtained is ^O.lmSl (FIG. 2). Note 
that this result does not crucially depend on the effec- 
tive statistics of vortices. We have also computed Rd 
by modeling the metallic phase as a classical hard-disk 
liquid of vortices 25 , and the resulting magnitude and the 
behavior of Rd is extremely close to the results we ob- 
tained within the vortex Fermi liquid framework 26 . 

At fields larger than the insulating phase or the MR 
peak, spinons (unpaired electrons) delocalize and im- 
pede vortex motion and suppress drag resistance. Fol- 
lowing Ref. 10, we use a semiclassical Drude formalism 
with statistical interactions between vortices, spinons, 
and Cooper pairs built in (e.g., the Magnus force on vor- 
tices when Cooper pairs move, etc.), and we find that 



with a finite spinon conductance, R s 1 > 0, 
Rd = Rd/ [(1 + Rv/R s ) 2 ] , 



(4) 



where R D is the vortex drag resistance with localized 
spinons, and R v = (h/2e) 2 a v is the vortex contribution 
to the resistance. Thus, when R s <C R v , the drag re- 
sistance is quickly suppressed to immeasurably small as 
spinon mobility increases. 

Lastly, we must estimate the vortex mass m v . Since 
there is still controversy over its theoretical value, we 
chose to estimate it from the experimentally measured 
activation gap in the insulating phase 5 ' 6 . When vor- 
tices condense, Cooper pair density fluctuations become 
gapped due to Higgs mechanism, and the gap which can 
be read off from the Lagrangian depends on m v . We 
conjecture that this gap is the activation gap. We find 



m v = 8Trh 2 n v T KT /E; 



gap J 



(5) 



with n v = B/(j) . For the InO film of Ref. 5 , T KT « 0.5K, 
and E gap « 1.6K at B = 9T, which implies m v w 
19m e , m e being the electron mass in vacuum. Note 
that Ref. 5,6 have reported the suppression of the ratio 
Tkt / E gap by increasing disorder. This is natural from 
(5), since Tkt I E gav ~ ^Tkt^v /n v , and disorder sup- 
presses both m v and Tkt- For comparison, this m v is 
close to some theoretical results m v ~ m e kpd for dirty 
superconductors 27 . All analysis above combines to yield 
the drag resistance behavior, which we plot in FIG. 2(a) 
for the film as in FIG. 2b of Ref. 5. 

Drag resistance within the percolation paradigm. 
Within the percolation picture, the non-monotonic MR 
arises from the film breaking down to SC and normal 
regions (described as localized electron glass) 12 . As the 
magnetic field increases, the SC region shrinks, and a 
percolation transition occurs. Once the normal regions 
percolate, electrons must try to enter a SC island in pairs, 
and therefore encounter a large Coulomb blockade ab- 
sent in normal puddles. The MR peak thus reflect the 
competition between electron transport though narrow 
normal regions, and the tunneling through SC islands. 
This picture is captured using a resistor network descrip- 
tion. Each site of the network has a probability p to be 
a normal (1 — p to be SC); each link is assigned a resis- 
tance from the three values Rnn, Rss, Rsn, that reflect 
whether the sites the link connects are normal (N), or 
superconducting (S) . An increase of the magnetic field is 
assumed to only cause p to increase. The important in- 
gredient is that Rs n has an activated form with a large 
gap representing the charging energy of the SC puddle; 
Rnn reflects the resistance of a localized electron glass 
with hopping conductivity, and Rss is mostly negligible. 

To calculate Rd, we first tune the parameters to make 
the single layer resistance resemble the experimental data 
in FIG. 2(b) of Ref. 5. Next, we place one such net- 
work (active layer) on top of another one (passive layer). 
Each link is treated as a subsystem, which might induce 
a drag voltage (an emf) e = IRd in the link under it in 
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the passive layer. When a link is between two normal 
(SC) sites, it is treated as a disorder localized electron 
glass (superconductor). The small resistance for the SC 
islands in this theory implies that vortices in the SC is- 
lands, if any, have very low mobility. We find that these 
vortices have negligible effect on i?c 28 , and Coulomb in- 
teraction provide the major drag effect (more theoretical 
details will be published elsewhere). Thus, two vertically 
aligned normal-normal (NN) links dominate the drag ef- 
fect. With the electron counterpart of (3) and the form 
of electron density response function from Refs. 29, we 
find Rd between two localized electron glass separated 
by vacuum is: 

7? ~ 1 R1R2 T 2 1 
D ~ 96tt 2 U/e 2 (e 2 nad) 2 2x [ ' 

Here, n « 5 x 10 20 cm~ 3 is the typical carrier density of 
InO 6 , d = 20nm is the film thickness, a — 25nm is the 
center-to-center layer separation, R\ y 2 are the resistances 
of the two NN links, x = a/ (2Tte 2 vd£, 2 L ) where v is the 
density of states and £ L w lnm is the localization length. 
In deriving Eq. (6) we used the averaged value of the 
inter-layer coulomb interaction along the z-direction of 
the layers. Solving the Kirchoff's equations for the two 
layers, we obtain the voltage drop and thereby the drag 
resistance. The results are shown in FIG. 2(b). 

Drag resistance in the phase glass theory. The phase 
glass model 13 focuses on the low field SC-metal transi- 
tion. It, therefore, does not allow yet a full calculation 
of the drag resistance. We leave a full analysis of drag 
within this theory for future work, but simply observe 
that in the glassy state the phases are ordered locally, 
and thus have no mobile vortices. The current-current 
coupling effects should therefore be absent, and the drag 
is mainly produced by Coulomb interaction. Thus we 
expect the sign of the drag voltage to be negative, and 
the drag resistance should be small due to the scarcity of 
excitations in a bosonic system. 



Summary. We have calculated drag resistance in bi- 
layer amorphous thin films separated by an insulator. 
Our calculation was carried out within the two compet- 
ing paradigms, vortex and percolation pictures, that may 
account for the phenomena observed at the breakdown 
of the superconductivity in amorphous thin films. In the 
percolation pictures, the drag resistance is due to inter- 
layer Coulomb drag and immeasurably small, ~ 10 _12 f2. 
In the vortex picture, however, the drag is caused by 
vortex motion. Since the vortex density is much lower 
than the charge carrier density, the drag resistance is 
orders of magnitude larger; our calculation shows that 
it reaches O.lmfi with the same sign as the single layer 
resistance. These estimates are made using parameter 
values that can easily be realized in experiments. Thus, 
the drag resistance measurement, albeit challenging due 
to the small scale of the maximum drag, can indeed pro- 
vide a sharp distinction between competing theoretical 
paradigms. In our future work, we will incorporate an 
interlayer Josephson coupling and analyze its effect on 
the drag resistance within the different paradigms. We 
expect that the drag resistance in both picture will be en- 
hanced, but the magnitude difference will remain. This 
would make the drag resistance easier to measure, and 
may not only improve its chances of determining the cor- 
rect theoretical paradigm, but also serve as a comple- 
mentary tool in the quantitative investigation of these 
fascinating systems. We note that yet another interest- 
ing possibility, which we leave for future research, is to 
enhance vortex drag by using high magnetic permeability 
insulators between the two layers. 
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